Probing the ^.^-wave Pomeranchuk instability by ultrasound 
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Selection rules of ultrasound attenuation and sound velocity renormalization are analyzed in view 
of their potential application to identify Pomeranchuk instabilities (electronic nematic phase). It is 
shown that the transverse sound attenuation along [110] direction is enhanced by the Fermi surface 
fluctuations near a d x 2_ y 2-wave Pomeranchuk instability, while the attenuation along [100] direction 
remains unaffected. Moreover the fluctuation regime above the instability is analyzed by means of a 
self-consistent renormalization scheme. The results could be applied directly to Sr3R,U2C>7 which is a 
potential candidate for a Pomeranchuk instability at its metamagnetic transition in strong magnetic 
fields. 
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I. INTRODUCTION 

In 1958, Pomeranchuk 1 considered stability conditions 
for normal isotropic Fermi liquids and argued that, be- 
sides a uniform (spin-dependent) deformation of the 
Fermi surface, i.e., an itinerant ferromagnetic instabil- 
ity, residual quasiparticle interactions can induce a non- 
uniform deformation of the Fermi surface. Since then, 
an inhomogeneous deformation of a Fermi surface caused 
by quasiparticle interactions is frequently referred to as 
"Pomeranchuk instability (PI)" It is sometimes called 
"nematic" electron liquid 2 because such a Fermi surface 
deformation can be characterized by a director order pa- 
rameter which is used to describe a nematic phase in 
conventional liquid crystals. The nematic order intro- 
duces an anisotropy in the otherwise isotropic physical 
quantities, and the strong anisotropy of the longitudinal 
resistivity observed in a clean two-dimensional electron 
system under high magnetic fields was interpreted as a 
signal of the development of the nematic order. 3 

In lattice systems the term PI is commonly used to 
refer to a Fermi surface deformation which breaks the 
symmetry of the underlying lattice. In case of two- 
dimensional square lattices, it was shown 4-7 that a PI 
with d x 2_ y 2 symmetry can be realized near van Hove fill- 
ings, which then reduces the C4 symmetry of the original 
Fermi surface to C2 symmetry. These findings have stim- 
ulated a number of further theoretical works 8-20 on the 
nature of a d x 2_ y 2-wave PI. Recently this phase has at- 
tracted a renewed interest in several fields of solid state 
physics. The occurrence of a d x 2_ y 2-w&ve PI is hotly 
debated as a possible explanation for the curious phase 
found in the bilayer ruthenate Sr3Ru2 07 under strong 
mag netic field. 21 ' 22 Also a dz,2_ y 2-wave PI may provide an 
explanation for the anisotropy in the magnetic excitation 
spectrum observed in cuprate superconductors, 23,24 and 
a similar argument is ongoing in the context of new iron- 
pnictide superconductors. 25,26 An analogous state in the 
spin channel was also proposed 27-29 and used to explain 
the so-called hidden-order phase found in URu 2 Si2 30 . 

Despite these extensive discussions, however, an iden- 



tification of such a d x 2_ y 2-w&ve PI is not an easy task. 
The problem lies in the difficulty in finding an experi- 
mental probe which couples to the order parameter of the 
d a; 2_j / 2-wave PI. Unlike a magnetic transition where we 
can detect the magnetic susceptibility using neutron scat- 
tering experiment, it is hard to detect the corresponding 
susceptibility for a d x 2_ y 2-w&ve PI which occurs in the 
charge degrees of freedom. In Ref. 31, it was suggested 
that a spatial pattern of the local density of states pro- 
vides a direct probe of the d x 2_ y 2-wave PI. It is however 
applicable only through surface probes, a scanning tun- 
neling microscopy, which is easily influenced by surface 
conditions. Therefore, it is desirable to have another way 
to identify the d x 2_ y 2-wave PI based on the bulk proper- 
ties of a sample. 

From the point of view that the characteristic 
anisotropy of a d x 2_ y 2-~wa,ve PI can be easily wiped out in 
a macroscopic scale by a formation of domains between 
two degenerate ground states [FIG. 1 (b) and (c)], 21,31 it 
had better concentrate on the fluctuation effects above 
the transition temperature 32,33 rather than concentrate 
on the phenomena deep inside of the ordered state. The 
basic idea presented in this paper is that we can use 
phonons for the detection of the d x 2 _ y 2 -wave PI. This 
idea is based on the fact that phonons not only have 
a coupling to the order parameter of a d I 2_ !( 2-wave PI 
through clectron-phonon interactions, but also have a 
propagation-direction-selective coupling. By combining 
this remarkable property of phonons with a well known 
fact that damping of phonons (sound attenuation) ex- 
hibits a remarkable divergence near a second order phase 
transitions, 34 we can pick up information on the d x 2_ y 2- 
wave PI. In other words, the fluctuation sound attenu- 
ation caused by a d x 2 _ y 2-wave PI can be regarded as a 
direct probe of the unconventional charge susceptibility 
characterizing this ordered state. 

In this paper, we investigate the behavior of the fluc- 
tuation enhancement of the sound attenuation caused by 
a dj.2_y2-wave PI. For this purpose, we use two differ- 
ent approaches. The first one is a phenomenological ar- 
gument using a Gaussian Ginzburg-Landau- Wilson ac- 
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FIG. 1: We consider a quasi two-dimensional system com- 
posed of a stack of two-dimensional square lattices. The or- 
der parameter of a d x 2_ y 2-wa,ve PI considered in this work is 
defined by Eq. (1) where r) x (rj y ) represents the strength of 
the bond along x (y) direction as shown in FIG. 1 (a). The 
two states in FIG. 1 (b) and (c) are degenerate where the 
thickness of arrows indicates the strength of the bond. 



tion. The second one is a microscopic analysis start- 
ing from a quasi-two-dimensional Hubbard model with 
on-site and nearest-neighbor repulsions, which is supple- 
mented by the self-consistent renormalization (SCR) the- 
ory 35 to treat the fluctuation effects. Both approaches 
are shown to give the same result that there is a selec- 
tion rule in the fluctuation enhancement of the sound 
attenuation and sound velocity softening near a d x 2_ y 2- 
wave PI, depending on the propagation directions and 
polarizations. This leads to a conclusion that we can uti- 
lize these properties to detect the d x 2__ y 2-w&ve PI with 
the aid of ultrasound measurements. 

The plan of this paper is as follows. In Sec. II, we 
present a phenomenological argument to see how the 
Fermi surface fluctuations near a d x 2_ y 2-wave PI man- 
ifest themselves in transverse sound attenuation. The 
main message of this paper is essentially given in this 
section. In Sec. Ill, we present a microscopic analysis 
of a rf x 2_ y 2-wave PI, and discuss the mean field phase 
diagram as well as the effects of fluctuations. Besides, 
we present a microscopic calculation of the fluctuation 
sound attenuation, which reinforces the result of the phe- 
nomenological argument in Sec. II. Finally in Sec. IV, we 
summarize our result and give a discussion on the appli- 
cation of our result to the bilayer ruthenate Sr 3 Ru207. 
We use the unit fcs = h = 1 throughout this paper. 



II. PHENOMENOLOGICAL APPROACH 

In this section, we present a phenomenological argu- 
ment to see how the Fermi surface fluctuations near 
a c? a .2_j ) 2-wave PI manifest themselves in phonon prop- 



erties. We pay special attention to the propagation- 
direction resolved sound attenuation. Because we con- 
sider the fluctuation effects above the transition point 
and mainly focus on the transverse sound attenuation, 
our argument is essentially different from that given in 
Ref. 10 where the longitudinal sound propagation below 
the transition point was discussed. 

The order parameter of a d^^-wave PI 4,5 which is 
considered in this work is defined by [see also Eq. (25)] 



V = \(Vx- Vy ) , 



(1) 



where f] x (rj y ) denote the bond strength along x (y) direc- 
tion (Fig.l) of the underlying square lattice. The order 
parameter n is odd under the permutation x <-> y. In the 
presence of nonzero r), the bond strength along x and y 
directions are different, hence it deforms the lattice and 
affects the phonon properties. 

Up to the lowest order in 77, the mode-coupling term 
<5Fgl between lattice (phonon) and the PI order param- 
eter is given by 

SF GL (U,n) = k [ U **(Q) ~ u vy(Q)]v(-QP) 

where Q — (q,Q z ) is a three-dimensional wavevector 
with the in-plane component q and the out-of-plane com- 
ponent Q z , k is the coupling constant, u xx — u yy = 
iq x Ux — i<ly u y with the lattice displacement field U = 
(u, U z ). The above form of the coupling is allowed from 
the symmetry reason because under the permutation 
x <-> y both n and u xx — u yy are odd, keeping the coupling 
invariant. 

We now study how phonons couple to the order pa- 
rameter of a d x 2_ y 2-wave PI depending on their propa- 
gation direction and polarization. For the moment, we 
consider the case where the wavevector and the polar- 
ization vector of a sound wave lie within a conducting 
layer, i.e., Q = (q, 0) and U = (u x ,u y ,0), because in- 
teresting results come out in this case. Then, the longi- 
tudinal (transverse) phonon mj, (ut) can be written as 
ul = q x u x + q y u y (ut = q y u x - q x u y ) with q = q/\q\. 

Consider a sound wave propagating along [100]- 
direction (i.e., q y = 0). In this case, using the fact that 
u xx — u yy = i\q\uL when q y = 0, we can write SFql in 
Eq. (2) as 



SF GL (u,n) = n^2(i\q\u L (q)jr](-q). 



(3) 



On the other hand, when a sound wave is propagating 
along [110]-direction (i.e., q x = q y ), the mode-coupling 
term (2) can be written as 



SF GL (u,n) = K^2(i\q\ur(qfjv(-Q)- 



(4) 



Equations (3) and (4) mean that, through the interaction 
(2), the longitudinal phonons couple to the d x 2_ y 2-w&ve 
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PI only when they have their wavevector along [100]- 
direction, while transverse phonons does only when they 
have their wavevector along [110]-direction. 

Using these results [Eqs. (3) and (4)], we next study 
the effect of the fluctuations of r) on sound attenuation. 
Dynamical behaviors of sounds are described by the fol- 
lowing action for phonons, 36 

SphO) = ^y 1 {y'L{q,u m )K L (q,\.u) m )u L {q,us m ) 



u* T (q, uj m )K T (q, iu m )u T (q, u r , 



(5) 



where the 
Ki°\q,iuj m 



has the form 



' : + (sl°') 2 q 2 with v being the 



unrenormalized kernel 

= "1 + (*i 0) ) 2 ? ! 
polarization index L or T. Here, s^' and the 
bare longitudinal and transverse sound velocities, 
Pion the mass density of ions, and u m — 2nTm is 
the bosonic Matsubara frequency. The information 
on phonon dynamics is extracted by studying the 
retarded quantity K^(q,u) = K v {q,\uo m — > u + i0+). 
In the presence of itinerant electrons, excitations of 
particle-hole pairs yield an imaginary part in of the 
form 38 lmK^(q,uj) = — u>/T p h, v (q) for small u). This 
results in a finite phonon lifetime T p h.v(q)- The sound 
attenuation a v {q) is related to the phonon lifetime as 
a v {q) = l/(s v T v h. iV ), or equivalently, 



a u (q) = —ImK^(q,w). 

StjU) 



(6) 



The phonon propagator [K^(q, oj)]^ 1 is modified in the 
presence of rj fluctuations. When we consider the Gaus- 
sian fluctuation region for simplicity, the corresponding 
action is given by 



where 



1 



r)(q,u> m ), (7) 



(8) 



is the d-wave density correlation function, yu cx T — T c0 
measures the distance from the mean field transition tem- 
perature T c o, £o and Td(q) — w §\q\ ar e the bare correla- 
tion length and the damping rate of fluctuations, respec- 
tively [for the microscopic definition, see Eqs. (41) and 
(42)]. 

After integrating out the rf x 2_ y 2-wave PI order param- 
eter r\ by performing the Gaussian integral in Sgl(v) + 
SS GL (u,ri) [where 6Sgl(u,t)) = J dr 5F GL (u,rj)}, we can 
show that the couplings [Eqs. (3) and (4)] give an addi- 
tional contribution 



4/Oion 



X d(q,iu m ) 



(9) 



to K v (q,'vjj m ). This leads to the divergent sound at- 
tenuation a v (q) oc — Imx^(<7, uj). Recalling that the in- 
plane longitudinal phonon (transverse phonon) couples 
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couple to r| 

divergent sound attenuation 
softening of sound velocity 

[100]: O 
couple to n 

divergent sound attenuation 
softening of sound velocity 

[110]: O 
couple to n 

divergent sound attenuation 
softening of sound velocity 

[100]: X 
not couple to n 

no divergent sound attenuation 
no change in sound velocity 



FIG. 2: Summary of the results derived in Sec. II. In the 
figure, ul(ut) represents longitudinal (transverse) phonons 
with the wavy line being the direction of the wavevectors. 



to i] only when the wavevector is along [100] direction 
([110] direction), we obtain 




«lioo oc ^ , 

Q!T110 OC /i^ 2 , 



(10) 
(11) 



while there are no divergent behaviors for sounds along 
the other directions, 



CXL110 

cmoo 



no divergence, 
no divergence. 



(12) 
(13) 



Here, L100 means the longitudinal sound wave propagat- 
ing along [100] direction, etc. 

In addition to Eq. (2) there is another relevant term, 

SF GL (u, V ) = K 'Yv*(q + q'Hq')^\q\uL(q)],(U) 

<3:<3' 

in which the longitudinal sound modes couple to the 
d x 2_ y 2-w&ve PI order parameter. As before, by integrat- 
ing out 7] in Sgl{v) + ^S' GL (u,r)) [where SS' GL (u,r]) = 
j dr 5F Glj (u, rj)] and assuming a three-dimensional be- 
havior of the fluctuations, we can show [see also the para- 
graph containing Eq. (63)] that this coupling gives the 
longitudinal sound attenuation the same divergent con- 
tribution as Eqs. (10) and (11), 



oll oc /j, ' 



(in all directions). 



(15) 



Note that this latter result is the same as that obtained in 
Ref. 39 for fluctuation sound attenuation near a magnetic 
transition, because the coupling (14) and the coupling 
discussed in Ref. 39 [Eq. (1) therein] have the same form. 

Before proceeding to the summary of this section, wc 
comment on the behavior of the sound velocity. The 



4 



anomaly in the sound attenuation is intimately related to 
the softening of the sound velocity. The renormalization 
of the sound velocity due to the fluctuation of rj is given 

by 

a v = 4 0) V 1 + (4 0) <?r 2 Re 6K?(q, uj^O). (16) 

This means that the divergence in the sound attenuation 
is accompanied by the reduction in the sound velocity 
(sound mode softening). 

The results obtained in this section is summarized in 
FIG. 2. Eqs. (10)-(13) and (15) mean that longitudinal 
sound attenuation always show a divergent behavior on 
approaching a second order phase transition of a d x 2_ y 2- 
wave PI. On the other hand, in-plane transverse sound 
attenuation are divergent only when their wavevectors 
are along [110]-direction. Likewise, longitudinal sound 
velocities always show a critical sound mode softening 
near a d x 2_j,2-wave PI, whereas in-plane transverse sound 
velocities show the softening only when their wavevec- 
tors are along [110]-direction. We therefore propose to 
detect the d a ,2_ y 2-wave PI by using a measurement of 
propagation-direction resolved transverse sound attenu- 
ation and sound velocities. 

III. MICROSCOPIC ANALYSIS 

In this section, we present a microscopic analysis of 
fluctuation sound attenuation near a d x 2_ J/ 2-wave PI in 
order to reinforce the phcnomenological argument given 
in the previous section. Firstly, we perform a mean field 
analysis of a d x 2_ y 2-wave PI starting from a single-band 
Hubbard Hamiltonian with on-site and nearest-neighbor 
repulsions, and draw the mean field phase diagram. Next, 
we investigate the fluctuation effect above the transition 
temperature by employing the self-consistent renormal- 
ization formalism. 35 Finally, based on the same micro- 
scopic model, we present a diagrammatic calculation of 
the fluctuation sound attenuation above the transition 
temperature, which is shown to give the same conclusion 
derived within the more phenomenological approach of 
the previous section. 

A. Mean field phase diagram 

We start from a single-band Hubbard Hamiltonian on 
a two-dimensional square lattice, 

H = Hq + H[ mp + Hi/ + Hy , (17) 

where 

H v = uY i n\{r i )n i {r i ) (18) 

and 

H v = v E EM^rvto) (is) 

(ri.rj) <r,a> 



describe the on-site and nearest-neighbor repulsions. 
Here n a (ri) = c' a (ri)c,j(ri) is the number operator for 
electrons on a lattice site Vi with spin projection a = ±1, 
and (ri,rj) means a bond between a nearest- neighbor 
lattice site r, and rj. 

In this work we consider the effect of an external mag- 
netic field _ff oxt perpendicular to the two-dimensional 
plane. Hence, the first term in Eq. (17) describes the 
kinetic energy plus the Zeeman energy, 

H = J^(e p - e F )n p , CT - fr^(n P;T - n Pi jJ, (20) 

P,t7 p 

where n P;(T = c p (7 c P;CT and c p , CT = -J= c CT (r t )c" ip - T '* 
with the number of lattice sites N. Here 
e p = —2t[cos(p x a) + cos(p y a)] — At' cos(p x a) cos(p y a) — 
2i"[cos(2p x a)+cos(2pj / a)] is the single-particle dispersion 
with hopping amplitudes —t, —t', —t" between nearest, 
next-nearest, and third nearest neighbors, eF is the Fermi 
energy under the zero magnetic field, and h = /iB-ffcxt 
with the Bohr magneton /Ub- Because we neglect the 
influence of magnetic field on orbital degrees of freedom 
(spin-orbit interaction, Landau diamagnctism) , the 
magnetic-field effect is absorbed into the spin-dependent 
Fermi energy €F t „(h) = ep + ah. Hereafter, Length and 
energy are measured in unit of the lattice spacing a and 
the nearest hopping amplitude t. 

The second term Hi mp in Eq. (17) represents the 
spin-independent short-range isotropic (s-wave) impurity 
scattering, 

ffimp = E Uim p( n ) n<7 ( ri )' ( 21 ) 

Ti , (7 

where the impurity potential M; mp obeys the Gaussian 

ensemble u imp (n) = 0, Ui mp (n)uimp(^) = "im P |u| 2 ^i,j 
with rtimp and u being the impurity concentration and the 
strength of the impurity potential. We treat the impurity 
potential using Born approximation, and the impurity- 
averaged bare Green's function is given by 

Gi 0) (p,ie„) = .„ * , (22) 

where £fc jCr = e p — cf.ct , and e„ = e n + T sign(e„) with the 
quasiparticle scattering rate T. 

As discussed in Ref. 7, Eq. (17) contains a phase with 
d a .2_j / 2-wave PI as a mean field solution. This phase ap- 
pears when the Fermi energy coincides with van Hove 
energy 5 eyH = 4(t' — t"). In this paper we consider a 
situation where the zero-magnetic-field Fermi energy eF 
does not satisfies the "van Hove condition" (ep 7^ £ vh), 
but a moderately large external magnetic field h tunes 
one of the spin-dependent Fermi energies ep.cr* (h) to a 
van Hove condition. In this situation, the spin compo- 
nent a* satisfying £f,o-*(^) = £vh is relevant to the oc- 
currence of the PI, and we hereafter consider only the 
"active" spin component a* . The antiferromagnetic cor- 
relations described by the on-site repulsion Hjj are sup- 
pressed because of the polarization under the strong ex- 
ternal magnetic field. Although the ferromagnetic state 
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FIG. 3: Schematic view of the density of states N(e) consid- 
ered in this paper. The density of states has a peak at van 
Hove energy evn = 4(t' — t"), and a spin-dependent Fermi 
energy e FiT touches eyH under magnetic fields. In the situa- 
tion shown here, the up-spin is assigned as the "active" spin 
component. 



could compete with the PI, we discuss the case where the 
PI is stabilized by choosing the parameters V ~ £7/2. 

We first make a mean field decoupling by taking the 
two contractions, 



H y = v E EVOMn)^)^-), (23) 

one by one. This yields the following mean field Hamil- 
tonian, 

H v F = -\Y. E m(ri)cUrMri + S)m 

Ti S=±x,y 

where rjs(ri) — V(cl, (ri+5)c a * (r^)}, and we have picked 
up only the "active" spin component a* satisfying the van 
Hove condition. We further assume that the parity of % 
with respect to 5 is even (i.e., r]s = rj-$). Then, we have 



rrMF 



\ E r ») + ^'( r *)) 4* (ri)Ca' (Vi + x) 
V(fi) + </(n)) c a* (r^Cj* (n + y) 



+H.c. 



(25) 



where r?(rj) = \[r) x {ri) - ?7 y (r,)] and r/ (r t ) = \[r) x {ri) + 
rj y (ri)]. Because rf possesses the four- fold symmetry of 
the underlying square lattice and can be regarded as a 
mass renormalization of the quasiparticlcs, we hereafter 
neglect it. On the other hand, a nonzero value of 77 cor- 
responds to the deformation of the Fermi surface which 
expands along the p x axis and shrinks along p y axis (or 
vice versa) and Eq.(25) describes a d x 2_ y 2-wa,ve PI. Go- 
ing to the momentum representation, we have 



rrMF 



(26) 



IP 



FIG. 4: The equation determining the impurity-averaged 
Green's function for the "active" spin component a*. The 
single solid line represent the impurity-averaged bare Green's 
function [Eq. (22)]. A double solid line, and dotted line with 
a cross represent the impurity-averaged Green's function, and 
impurity scattering, respectively. 



Several comments arc in order. In a situation where 
the PI could be stabilized under the zero magnetic field, 
it would compete 40 with the so-called d-density wave 
state, a state with rj(q = ±(tt,tt)). Under a moder- 
ately large magnetic field, however, the d-density wave 
state is not expected to be competitive anymore since 
this state shows diamagnetic properties 41 suggesting that 
such a state is suppressed under sizable magnetic fields. 
Solving these problems is beyond our scope because it 
would require to include Landau diamagnetism such that 
the resulting analysis becomes considerably more compli- 
cated. Therefore in the following we assume that a uni- 
form rj(q = 0) gives the ground state, and the momentum 
q appearing in Eq. (26) is small (i.e., q <C 1). 

The equation determining the impurity-averaged 
Green's function for the "active" spin component a* is 
diagrammatically shown in FIG. 4 where the single solid 
line represents the impurity-averaged bare Green's func- 
tion [Eq. (22)]. Hereafter, we neglect the last diagram 
in FIG. 4 assuming that the impurity potential is suffi- 
ciently weak. Due to the symmetry difference between 77 
and the impurity potential, this approximation is exact 
up to the linear order in r\. Physically this means that the 
impurity effect on the second order transition line is rig- 
orously treated while the first order transition line may 
slightly deviate from the exact result. We note here that 
even with this simplification, the required full numerical 
calculation is rather involved. With this approximation, 
the equation for the impurity-averaged Green's function 
is solved as 



(27) 



where the order parameter r] is determined by the follow- 
ing self-consistent equation 



7] = V d p f r {£p,a* -d p r]). 
J P 



(28) 



Here we have introduced a shorthand notation f 

j p 

~k Sp> an d the quantity 



where d p — cosp x — cosp y , and f](ri) = J2 q <7(< / )e 1<J "" i . 



(29) 
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may be considered as the impurity-averaged Fermi dis- 
tribution function. The corresponding mean field free 
energy is given by 



N 



V 

2V 



T E/ ln 



\G^(P^e n )\ 2 



(30) 



In the pure limit (r — > 0), the function /r be- 
comes equal to the Fermi distribution function /o(a;) = 
(1 + e ^/ T )- 1 due to the identity T^ £ „ G ( °){p,ie n ) = 
/o(Cp,<r*)- This would simplify the self-consistent equa- 
tion (28), and the free energy (30) would be re- 
duced to a simpler form ^ = ^ - T / In | (1 + 
e -(£p-d P r))/Ty^ _|_ e -W T )|. However because we con- 
sider the effect of non-magnetic impurities in this work, 
we need to use Eqs. (28) and (30). 

The mean field phase diagram can be drawn by making 
a Landau expansion 17 of the free energy (30) as 



N 



N 



V 2 ) 



\ 4 



where the coefficients ctQL and 6ql are given by 



&gl 



^-j4{-Mz P ,*-j), 
| jf 4 (/?'(&.»•)). 



(31) 

(32) 
(33) 



with a convention /f^p.a) = ~q^~^ fv (£p,<r)i etc. 

A condition ogl = defines a second order transition 
line T c o(h) in the mean field approximation, if the coef- 
ficient of the quartic term of the GL free energy (&gl) is 
positive. In case of a negative quartic term the transition 
becomes of the first order, and the first order transition 
line T*(h) is determined by the conditions Fmf(v) = 
and dFyiF^/dri = for a nonzero r\ where Fmf is de- 
fined in Eq. (30). 

We calculate the mean field phase diagram by using a 
square mesh of 500 x 500 in the Brillouin zone. In FIG. 5 
(a), we show the calculated mean field phase diagram for 
a moderately clean sample (r = 0.01). Not surprisingly, 
the phase diagram in this clean case is quite similar to 
that already obtained in the pure limit (r = 0) in the pre- 
vious work, 42 since as will be discussed below Eq. (37) 
in the next section, our model is similar to that used in 
Rcf. 42 except for the difference in a spin mixing. The 
important feature in the phase diagram is that the tran- 
sition is of second order at higher temperatures while at 
lower temperatures a first order transition occurs. In the 
figure, we used a non-zero third-nearest neighbor hopping 
t" = 0.2. This is because when we use smaller values of 
t", a finite momentum state rj(q ^ 0) is stabilized due to 
a nesting condition. 

Next let us see how a slight increase in the impurity 
scattering modifies the mean field phase diagram. 43 In 
FIG. 5 (b), the mean field phase diagram for a dirtier 
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FIG. 5: Phase diagrams obtained from Eq. (17) for two dif- 
ferent strengths of the impurity scattering, T = 0.01 [(a)] and 
r = 0.06 [(b)]. Note that length and energy are measured in 
unit of the lattice spacing a and the nearest hopping ampli- 
tude t. We use U = 3.0 and V = 1.5 for electron repulsions, 
and t' — 0.4 and t" = 0.2 for hopping. Two (thick and thin) 
dashed lines are defined by ogl = where the thick dashed 
line represents a mean field second order phase transition line 
T c o(h) while the thin dashed line is an unphysical one be- 
cause it is replaced by the following first order transition. A 
line with full circles show a first order transition line T*(h) 
defined by Fmf(ti) = and (9Jmf(t/) / dij = 0. A line with 
open circles represents the fluctuation renormalized second 
order transition line T c (h). A dotted line represents a tem- 
perature below which the quartic term &gl [Eq. (33)] of the 
Ginzburg-Landau free energy becomes negative. 



case (r = 0.06) is shown. Firstly, we see that the PI is 
easily suppressed by the increased impurity scattering. 
This is understandable because the order parameter of 
the PI and the impurity potential have different sym- 
metry. Hence, the underlying mechanism is analogous 
to what happens to a disordered d-wave superconductor 
where Anderson's theorem is violated due to the uncon- 
ventional nature of the order parameter. Secondly, the 
first order transition line at low temperatures disappears 
while the second order transition still survives. These 
two features concerning impurity effects are important 
when we discuss the relevance of the Pi-scenario to the 
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metamagnetic transition found in Sr 3 Ru20y in Sec. IV. 



B. Self-consistent Renormalization (SCR) 
treatment of the Fermi surface fluctuation 

In this subsection, we study how the mean field result 
obtained in the previous subsection is modified by fluc- 
tuation effects. Following Refs. 14 and 18, we define the 
(thermal) d-wave density correlation function 



1 

N 



l/T 



where n d (q) = J2 P d P c P+q / 2 ,cr* c p-q/2,<?' is tnc d-wavc 
density operator in the "active" spin component a*. In 
general, \d has the following structure 



1/Xd 



(35) 



where xT i s the irreducible part of Xd- To proceed fur- 
ther to the detailed calculation of Xd, it is convenient to 
rewrite the interaction Hamiltonian Hy in the momen- 
tum space, 



H 



y - 2^ E E^- 



P2 C p! ,(X C P2 ,CT 



Pl>P2,P3 CT,Cr' 



XC p3,<r' C Pi+P3-P2,a-'j (36) 



as Vp 1 — p 2 
where d p 



where Vp = 2V(cosp :E -|-cosp, / ). Now we decouple V Pl - P2 

V[d pi d p2 + s pi s p2 -\-ppiPp-} + Ppi^Pp2^]> 
cosk x — cos k y , s p = cosk x + cosk y , and 

Pp^ = sinfca; ± sinfc y . Because the four van Hove points 
in the Brillouin zone enhance the interaction with d x 2_ y 2 
form factor 5 d p , we hereafter consider only this channel. 
After changing the ordering of fcrmion operators and set- 
ting pi = p + q/2,p 2 = p' + q/2,p 3 = p' -q/2, we have 



V 



4 = -^E^(^(-«)' ( 3? ) 



This is essentially the same model as used in Refs. 12 and 
42 except for the fact that the interaction H' v acts only 
among the active spin component a*. 

The so-called random phase approximation (RPA) for 
Xd is obtained when we adopt the simplest building block 



x" (g,iw m ) 



E / 



i£n) 



xG ( °)(p + q,ie n + iLu m ) (38) 



as the irreducible x 1 " ■ Hereafter, we write the <i-wave 
density correlation function in a dimcnsionless form by 
introducing Xd = Xd/xTo- Then, the RPA d-wave den- 
sity correlation function is given by 

l/Xd,RPA(q,iw m ) = l-Vxd*o(q,iUm)- (39) 



+ 



(c)(.irr 

Xd 




FIG. 6: Diagrams for the fluctuation propagator of d x 2_ y 2- 
wave PI [(a) and (b)] and the irreducible d x 2_ y 2-w&ve correla- 
tion function giving the mode coupling between fluctuations 
[(c)]. The solid line represents the impurity-averaged bare 
Green's function of quasiparticles, the dashed line the next 
nearest neighbor repulsion [Eq. (36)]. 



Near the second order PI, the retarded function 
X^rpa(9' w ) = Xd,RPA(q,iu m -> w + iO+) takes the same 
form as Eq. (8), 



w {q)\q\ 



(40) 



where fi = Vo-gl oc T — T c0 with a GL and T c0 being 
the quadratic coefficient in Eq. (32) and the mean field 
transition temperature. Also, £q and wo are given by (see 
Appendix A) 



e 



V d. 



r " T{ ^'\vl+vl) 



+ 



12 

S " t{ ^'\ m x + M v ) , (41) 



wo(q) 



= V d; 



/f(^p,c*) 



p p \2T/q c + iv Pta * q) ' 



(42) 



where Mj = (l/2)d 2 ^ p /dkj for j = x,y and q c is the 
wavevector cutoff. We note that, as was already men- 
tioned above Eq. (27), neglecting the impurity vertex 
corrections, 44 which would otherwise transform the dy- 
namical behavior into a diffusive one, remains exact in 
the calculation of these quantities. 

In the presence of interlayer coupling, the critical be- 
havior is expected to be governed by three-dimensional 
fluctuations, 45 and we assume in this work an anisotropic 
three-dimensional (3D) behavior of the fluctuations. Mi- 
croscopically, the three dimensionality is introduced by 
adding a z-axis dispersion 5e Pz = —5tcos(P z ) into 
Eq. (20). Instead of doing so, however, we here intro- 
duce the three dimensionality in a more phenomeno- 
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logical way, by replacing the two-dimensional wavevec- 
tor q appearing in XdiQi'^m) with a three-dimensional 
wavevector Q' = (q,Q z /l), namely, 



Xd(Q',iu m ), 



(43) 



where 7 = £o/£cO is the anisotropy parameter given by 
the ratio of the in-plane correlation length to the out-of- 
plane correlation length. The above procedure is justi- 
fied at least for a system with anisotropic ellipsoidal 3D 
Fermi surface. Here we would like to mention the recent 
proposal in the context of high-T c cuprates 46 that the 
intcrlayer-configuration of the nematic order becomes al- 
ternate pattern of 77 > [FIG. 1 (b)] and r\ < [FIG. 1 
(c)]. We do not take into account such a possibility here, 
since this configuration costs the elastic energy of the 
crystal lattice. We would like to point out that, while a 
minor change in the form of the propagator in Eq. (40) 
is needed if the cL^^-wave PI has such a configuration, 
the main result of this paper is not changed, at least 
qualitatively. 

Now we consider the effect of mode coupling between 
the fluctuations. In the bosonic languages, 47 ' 48 this is 
given by the self-energy renormalization for the fluctu- 
ation propagator coming from the quartic term of the 
action. In the fermionic languages which we adopt in 
this work, the corresponding mode coupling is given by 
the three diagrams 49 shown in Fig. 6, which give 



x7 = x%+s^\ 

5x7 = -3b GL V 2 TY / 



jd 3 Q' 
(2tt) 3 



(44) 

Xd(Q',iw m ),(45) 



where &gl is the coefficient of the quartic term in the 
Ginzburg-Landau free energy given by Eq. (31). Substi- 
tuting Eq. (44) into Eq. (35) and expanding it in terms 
of 5x7 ■> we obtain 



1/Xd = VXd,RPA - 5x7 



(46) 



where 5x7 a ^ so depends on Xd through Eq. (45). 

As discussed in Rcf. 48, the SCR formalism is a self- 
consistent one-loop approximation for the mass renor- 
malization 



/1 = i/xd(Q' = o, 



0). 



(47) 



Substituting Eq. (47) into Eq. (46) and making the con- 
tour deformation for the Matsubara summation, we ob- 
tain the following SCR equation to determine /x, 



, ,r* ^ f 



(48) 



(27T) 

xImx£(Q», 



\2T) 



(49) 



where Q' c is the wavevector cutoff. The frequency in- 
tegral in Eq. (49) can be performed using the relation 




0.06 0.08 0.10 

T Tco=0.0776 



FIG. 7: Temperature dependences of parameters character- 
izing the bare d-wave density correlation function; the mass 
term /io, correlation length fo, and damping coefficient wo. 
The mode-coupling parameter &gl is also plotted. The pa- 
rameters used in this figure are the same as in FIG. 5 (a). 



coth(a;/2) = l + 2/(c z -l) and dx(e x - l)" 1 [a;/(a; 2 + 
c 2 )] = (l/2)[ln(c/27r) - tt/c - *(c/2tt)]. Then, introduc- 
ing the normalized wavevector Q' = £,oQ' , the final result 
can be written as 



S/j, 



x2Z~( \nZ~ - 

Q V Q 2Z~ 



where *&(Z) is the digamma function, /z~ = ji + Q /2 , 
Z~, = fj,~ w \Q / \/(2irT( ). When the condition Z~ < 
1 is satisfied, we can set 2Z~,[lnZ~, — (1/2)Z~, — 
^(Z~, )] — > 1 and the above expression reproduces the 
result of the Hartree approximation for the mass renor- 
malization due to classical fluctuations, 37 which means 
that Eq. (50) can describe both quantum and classical 
fluctuation regions. 

In FIG. 7, we plot the temperature dependences of 
several parameters characterizing the bare d-wave den- 
sity correlation function, the mass term /x , correla- 
tion length £ ) an d damping coefficient w , as well as 
the mode-coupling parameter 6gl- Strictly speaking, 
wo(q) in Eq. (42) depends on q but we average them 
using two values at q = x and q = (x + y)/2, i.e., 
wo(q) -» < v = (l/2){iuo(x)+iuo[(8 + y)/2]}. From the 
figure, we can see the following three points. First, as was 
already discussed in the pure limit 42 (r = 0), the bare 
mass term fiQ is as small as <~ 0.1 in a wide temperature 
region above the mean field transition temperature T c0 . 
Second, the correlation length is shorter than the lattice 
spacing a = 1, which can be understood from Eqs. (41) 
because the quasiparticle velocity |u p | is quite small near 
van Hove filling. Third, the mode coupling parameter 
6ql is one order of magnitude larger than other param- 
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0.05 0.07 0.09 0.11 0.02 0.03 0.04 0.05 

T T 



FIG. 8: Temperature dependence of the renormalized mass 
fi — l/Xd{Q' = 0,iuj m = 0) (open circles) and the bare 
mass fxo = l/Xd,RPA(Q' = 0,io> m = 0) (dashed line) for 
T = 0.01 [(a)] and T = 0.06 [(b)]. We use the wavevector cut- 
off Q' c — 0.0555 and anisotropy parameter 7 = 5, and other 
parameters are the same as in FIG. 5. Two curves are calcu- 
lated along a fixed magnetic field h = hyu where the mean 
field transition is always of second order. The fluctuation- 
renormalized transition temperature T c is determined by a 
condition fi(T c ) < 10~ 3 . 



eters (po,S,o, V^o) characterizing Gaussian fluctuations. 
All these features result in strong fluctuations near the 

In FIG. 8, we plot the temperature dependence of the 
renormalized mass /i = l/xd(Q' = 0,iuj m = 0) as well 
as the bare mass po — 1/X(2,r.pa(<3' = 0, \uj m = 0). As 
expected, p{T) coincides with po(T) well above the mean 
field transition temperature T c q but it deviates from the 
mean field result near and below T c0 . From the figure we 
determine the fluctuation-renormalized transition tem- 
perature T c by a condition p{T c ) < 10 -3 . In FIG. 5 we 
have also plotted T c {h) thus determined. We can see that 
the fluctuation region becomes smaller on approaching 
the first order transition line. 



C. Fluctuation sound attenuation 

We consider the clcctron-phonon interaction following 
the argument given by Walker, Smith and Samokhin, 50 
and begin with the tight-binding Hamiltonian 

ffd-ph = ^2g{5)^2(u{ri) -u(ri + 5)j 

x (clin^in + 5) +4(r i + 5)c CT (r i )) (51) 

where u(ri) is the displacement of the ion at r^, 8 is the 
lattice vector, g(S) = X7t(5) with the hopping amplitude 
t(S) between site rj and Vi+6. Following the procedure in 
Ref. 50 and neglecting umklapp processes, the interaction 
between electrons and a sound wave with wavevector K 



and polarization v is given by 

I JO) 

^r B K:" C P + K,* C P;°, ( 52 ) 

where B K u — b K v — b*_ K v with the phonon annihila- 
tion operator bjc,v, &v is the polarization vector, and the 
electron-phonon vertex function F^g ^ (P) is given by 

xcos(P-S). (53) 

Here, all momenta denoted by capital letters are 3D 
vectors. However, we consider below the case where 
the wavevector and the polarization vector of a sound 
wave lie within a conducting layer because interesting 
results come out in this case, and we set in the fol- 
lowing P = (p, 0), K = (fc,0), etc. In this work 
we consider the clcctron-phonon coupling through at 
most the next nearest neighbor interactions. Assum- 
ing that g(S) is proportional to <5, we make an expan- 
sion g(S) = g n .n.8 n .n. +ffn.n.n.^n.n.n. where 

yn.n. lyn.n.n. 

the coupling constant through the nearest (next nearest) 
neighbor interactions. 

Then, for transverse phonons propagating along the 
[100] direction, because the coupling via the nearest 
neighbor interaction disappears in this case, the domi- 
nant coupling is given by the next nearest neighbor in- 
teractions as 

Ftioo(p) = T T ioo(2 sin p x sin p y ), (54) 

while for transverse phonons propagating along the [110] 
direction, 

Ftuo(p) = T T110 {cosp x - cosp y ), (55) 

where T T i 00 = -gn. n.n.\/L z /p ion (sTwo) 2 [T T no = 
—gn.n\/ L z / pi on (sTiw) 2 ] with system's c-axis dimension 
L z . The electron-phonon vertex function for the longitu- 
dinal phonons is given by 

Fl(p) = T L {q 2 x cos p x + q 2 y cos p y ) (56) 

with Tl = —g n .n. \/L z / Pion(sL) 2 - The important point 
here is that in each case F~ v }{p) can be expressed as 
^{K v}(P) = TW P with a form factor W p and the cor- 
responding electron-phonon coupling constants T. For 
T100 phonons, W p is equal to the dry-symmetry form 
factor d' p = 2 smp x sinp y , while it is equal to the d x 2_ y 2- 
wave form factor d p for T110 phonons. In case of lon- 
gitudinal sounds, W p = s p + dp for L100 phonons while 
W p = s p for LI 10 phonons where s p is defined below 



Eq. (36). We note that the inclusion of an isotropic vertex 
function into Eq. (56), which ensures the charge neutral 
condition 50 for longitudinal phonons, does not change the 
main result discussed below. 

Using the electron-phonon interaction (52), the sound 
attenuation is calculated from the phonon Green's func- 
tion 

[D„(fc,ifim)] _1 - [Dp^k^Slrn)]- 1 - ILyfciSlrrJpr) 

where D?\k,iQ m ) = -(4°>)7Rl + (4°l) 2 ] * the 
bare phonon Green's function, and n„(fe, ifl m ) is the 
phonon self-energy. The retarded phonon self-energy, 
n^(fc,i£! m ) = n„(fc,if2 m — > fl + i0+), determines the 
phonon lifetime as l/r p h,i/ = — (cj^I) 2 j^lmllfj-(k, fi), and 
the sound attenuation a„(fc) for the phonon with wave 
vector k and polarization v is given by 

f(J 0) ) 2 \ 1 

M*) = - ( -jgf- j ^imn*(fc,n). (58) 

Strictly speaking, the sound velocity in Eq. (58) should 
be understood as a renormalized one, 

s v = sW^/l-Ren*(fc,n-»0), (59) 

with the bare velocity si°\ However, because the renor- 
malization of the sound velocity due to the collective fluc- 
tuation is usually quite small in the experimental res- 
olution 56 [{s v - si 0) )/ s ^ 0) ~ 10_3 - 10" 4 ], wc use the 
bare sound velocity in Eq. (58). Note however that the 
unrenormalized sound velocity should be used only in 
Eq. (58). As mentioned at the end of Sec. II, the mea- 
surement of the softening of the sound velocity is also a 
key experiment to detect the d^^-wave PI. 

Consider first the sound attenuation caused by itiner- 
ant electrons in the absence of the d x 2_ y 2-wave PI. In 
this case the irreducible phonon self-energy is given by 38 
[FIG. 9(a)] 

n(°)(fe,iO ro ) = T 2 T]T / W 2 G^(p,ie n ) 

This gives the sound attenuation without 
the Fermi surface fluctuations, a]f F (k) = 

-[(4°^)V^imnW*(fc,ft). 

Now we consider the effect of the d^^-wave PI. The 
relevant diagrams are shown in FIG. 9(b)-(f) where we 
have neglected impurity vertex corrections due to the rea- 
son mentioned above Eq. (27). The simplest diagram 
is given by FIG. 9(b) which corresponds to the sound 
attenuation derived from the mode-coupling term (2). 
Diagram (b) gives the sound attenuation derived from 
the mode-coupling term (14), and is analogous to the 
Aslamazov-Larkin diagram 51 for the fluctuation conduc- 
tivity in a superconductor. Here we have neglected the 
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FIG. 9: Diagrams for the phonon self-energy ID(fc, if2 m ). An 
open circle represents the electron-phonon vertex in Eq. (53), 
a full circle the d a .2_ !/ 2-wave vertex d p , a wavy line the phonon 
Green's function, and a double zigzag line the fluctuation 
propagator defined in FIG. 6. 

vertex corrections due to mode-coupling [e.g., FIG. 1(c) 
in Ref. 52]. Diagrams (d)-(f) give the contributions which 
are not described by the phenomenological (bosonic) ar- 
gument in Sec. II, and analogous to the Maki-Thompson 
diagram 53 ' 54 [diagram (d)] and density of states dia- 
grams 55 [diagram (e) and (f)]. Note however that the 
analytical expression of each diagram is different from 
that of a superconducting fluctuation contribution be- 
cause the PI is a kind of diagonal long range order while 
the superconductivity is an off-diagonal long range order. 

For the moment, we focus on the transverse sound at- 
tenuation because it is this case where the fluctuation 
sound attenuation is propagation-direction selective. We 
first consider diagram (b). For the transverse sound 
along [100] direction, because the electron-phonon ver- 
tex function Ftioo(p) [Eq. (54)] has the (^-symmetry, 
the phonon self-energy does not couple to the collective 
fluctuation of the d a .2_ y 2-wave PI. Hence the phonon self- 
energy has no divergent contribution, and is given by 

n T io (fe,ifi m ) = U^l Q0 (k,iQ, m ), (61) 

where n^ 00 is defined in Eq. (60) with W p = d' p . On 
the other hand, the transverse sound along [110] direction 
does couple to the d x 2_ y 2-wa,ve Fermi surface fluctuation 
because the electron-phonon vertex function Ftho(p) is 
proportional to the c? :r 2_ J) 2-wave form factor. This gives 
the divergent phonon self-energy 

n Tno (fe,iO m ) = -(T T ii ) 2 Xd(fe,i^m), (62) 

where Xd{k, ifi m ) behaves as Eq. (8) near the second or- 
der transition. 
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(T-Tc)/Tc 






(T-Tc)/Tc 



0.8 

(T-Tc)/Tc 



FIG. 10: Calculated ultrasound attenuation for T100 phonons 
(open circles) and T110 phonons (full circles) for T = 0.01 
[(a)] and T = 0.06 [(b)] at constant field h — hvH- The param- 
eters used in this figure are the same as in FIG. 5 and FIG. 8. 
The data are normalized by a MF at the mean field transition 
temperature T c q where q mf is defined below Eq. (60). 



Next we consider diagram (c). In this case, it is easy 
to show that the triangle blocks in this diagram van- 
ish for both T100 and T110 phonons, hence there is no 
contribution from this diagram. This is consistent with 
the argument in Sec. II that Eq. (14) couples only to the 
longitudinal sound, because diagram (b) gives the sound 
attenuation derived from this mode-coupling term (14). 

We then consider diagrams (d)-(f) which were not 
taken into account in the phenomenological argument in 
Sec. II. Note that if we pick up only the contribution 
just at the Fermi surface assuming a constant density of 
states, then we can show that diagram (d) and diagrams 
(c) and (f) cancel. However, once we include the energy 
dependence of the density of states, there is a nonzero 
contribution from these diagrams as in the case of the 
quartic term of the GL functional. By power counting, 
we can show that these diagrams give a logarithmically 
divergent sound attenuation olt ~ In// where // is the 
renormalizcd mass term (47) . We think that in the actual 
experiment this logarithmic divergence is negligible com- 
pared with the stronger divergence caused by Eq. (62). 



FIG. 11: Calculated ultrasound velocity for T100 phonons 
(open circles) and T110 phonons (full circles) for F — 0.01 
[(a)] and T = 0.06 [(b)] at constant field h = h VH - The 
parameters used in this figure are the same as in FIG. 5 and 
FIG. 8. The data are normalized by the bare sound velocity 
s(°\ Here, we have used a moderately large electron-phonon 
coupling constant T 2 = 0.001 to enlarge the signal, whereas 
the actual softening would be much smaller. 56 



Now we come to the case with longitudinal sound 
modes. Concerning the fluctuation contribution to the 
longitudinal sound attenuation, it is sufficient to con- 
sider diagram (c) in FIG. 9 irrespective of the propa- 
gation direction because as you can see in the following 
this diagram always gives the divergent behavior with the 
same exponent. After analytical continuation, diagram 
(c) gives 



n?(fc,fi) 



-iQB 



dw 



sinh 2 (^y) 



rf 3 Q'[lm X J(Q',^3) 



where B = r L J2 E J p W p d 2 p [G^(p,ie n )} 3 = 
~2~ J p Wpdpfr(£p,a*) comes from the triangle block 
in the diagram and Q' = (q,Q z /j) is the rescaled 
three-dimensional wavevector introduced in Eq. (43). 
Roughly speaking, B is proportional to the degree 
of particle-hole asymmetry, and in case of a perfect 
particle-hole symmetry this contribution vanishes. By 
power counting, we can show that this gives a divergent 
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sound attenuation ~ \i~ 2 where /i is the renormalized 
mass term (47). Therefore, as was already discussed in 
Sec. II, the longitudinal sound attenuation always show 
a divergent behavior on approaching a second order PI. 

To summarize, a microscopic calculation of the fluctua- 
tion sound attenuation in this subsection gives the same 
result as obtained by a phenomenological argument in 
Sec. II. Results of this subsection coincide with the main 
results in Sec. II, i.e., Eqs. (11), (13) and (15), if we re- 
place the bare mass (j,o with the renormalized mass \x. 

In FIG. 10 we plot the transverse sound attenuation 
along [100] and [110] directions as functions of temper- 
ature. We also plot in FIG. 11 the transverse sound ve- 
locities along [100] and [110] directions as functions of 
temperature. As was already discussed in Sec. II, there 
is a fluctuation contribution to [110] phonons while there 
is no contribution to [110] phonons. Comparing data for 
a cleaner system [FIG. 10(a)] to that for a dirtier sys- 
tem [FIG. 10(b)], we see that the dirtier system shows 
broader fluctuation behavior in the reduced temperature 
(T - T c )/T c . 



IV. DISCUSSION AND CONCLUSION 

We have studied the effect on sound properties of the 
Fermi surface fluctuations near a 2 -wave PI, and 

discussed that there is a propagation-direction-dependent 
selection rule in the fluctuation transverse sound attenua- 
tion and sound velocity softening. As is shown in FIG. 2, 
FIG. 10, and FIG. 11, the transverse sound attenuation 
and sound velocity softening along [110] direction arc 
enhanced by the Fermi surface fluctuations while those 
along [100] direction are not affected. Also it was ar- 
gued that there are always fluctuation contributions to 
the longitudinal sound attenuation and sound velocity 
softening. We note that a qualitatively similar conclu- 
sion can be reached for a second order structural transi- 
tion which breaks the lattice symmetry in the same way, 
but is not caused by a Fermi surface instability. Such a 
transition would not be identified as a genuine Pomer- 
anchuk instability. 

As was already mentioned in Sec. I, the possibility of 
the d z ,2_ y 2-wave PI discussed in this paper has been de- 
bated 16,21, 22,31 ' 42 ' 57,58 as a possible explanation for the 
anomalous phase found in Sr3Ru207 under strong mag- 
netic fields, and it provides us with a good opportunity 
to apply our results. In earlier experiments for this mate- 
rial, a metamagnetic transition was found 59 ' 60 at a mag- 
netic field B m where the magnetization shows a steep 
jump. Since several non-Fermi-liquid properties are ob- 
served around B m , this system was originally discussed 
in terms of the metamagnetic quantum critical point, 61 
and the importance of van Hove singularity for the meta- 
magnetic transition was discussed. 62 However, the sub- 
sequent experiments using an ultra-pure sample found at 
least two transitions, 63,64 and later on it was revealed 21 
that the two transitions reflect the boundary of a new 



ordered phase which is accompanied by metamagnetic 
transitions. Based on a consideration on resistivity data, 
it was proposed 21 and demonstrated 16 that the observed 
new phase can be an ordered state with d a .2_ y 2-wave PI. 

On this background, it is tempting to compare 
our results with the experiments for bilayer ruthenatc 
Sr 3 Ru207. Firstly, as was already discussed in Ref. 42, 
the calculated phase diagram [FIG. 5 (a)] looks quite 
similar to what was observed in experiments (Fig. 3 in 
Ref. 21). Secondly, the earlier experiments can be inter- 
preted as a result of the impurity effects which narrow 
the area of the ordered phase [FIG. 5 (b)], because two 
closely located transition lines (beyond the experimen- 
tal resolution) would look like a single transition line in 
experiments. Thirdly, the non-Fermi-liquid behaviors ob- 
served in resistivity, 59,60 ' 65 specific heat, 60,65 ' 66 and ther- 
mal expansion 67 can be interpreted as a result of the 
Fermi surface fluctuations near the d x i ^ y 2 -wave PI dis- 
cussed in Sec. Ill B. This argument can also be applied to 
the non-Fermi-liquid behavior observed in nuclear relax- 
ation rate, 68 if we take into account the spin-orbit scat- 
tering. These facts suggest that the d a .2_ y 2-wave PI is a 
promising scenario for the curious phase observed in this 
material. Hence, it is interesting to apply our main re- 
sult to Sr 3 Ru 2 7 , i.e., Eqs. (11), (13) and (15) [note that 
Eq. (12) is replaced by Eq. (15)], and we propose to mea- 
sure the propagation-direction resolved transverse sound 
attenuation and sound velocity softening in Sr3Ru207, 
since it can determine experimentally the presence or ab- 
sence of the PI in this material. 

Finally we comment on the recent papers 69,70 on 
the microscopic mechanism of the ei^^-wave PI for 
Sr 3 Ru 2 07. In these articles, the authors point out the 
importance of the quasi-onc-dimcnsional ruthenium or- 
bitals (d xz and d yz ) for the occurrence of the d x 2_ y 2- 
wave PI. On the other hand, our microscopic analysis 
in Sec. Ill essentially uses quasi-two-dimensional ruthe- 
nium orbital d xy . We point out that while the micro- 
scopic analysis in Sec. Ill needs small modifications 71 if 
we start from d xz and d yz orbitals, the main conclusion 
remains unchanged that the transverse sound attenuation 
along [110] direction and the longitudinal sound attenu- 
ation in all directions are enhanced by the Fermi surface 
fluctuations while the transverse attenuation along [100] 
direction is not affected. This can be inferred from the 
phenomenological nature of the argument given in Sec. II 
to arrive at Eqs. (11), (13) and (15). 
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APPENDIX A: CALCULATION OF ^ R (q,ui) 

To evaluate x^ RPA (q,w), we need to calcu- 
late xYQ R {q,uj). On Matsubara (imaginary) axis, 
Xdo(9) it^m) is given by Eq. (38). We transform the Mat- 
subara sum into an integral over the contour as shown 
in FIG. 12 using T£ £ ^(ie„) = § ^ tanh(^)F( £ ), 
and perform analytical continuation iuj m — > u) + i0 + . 
Contour C\ and C3 give the real part of x^ (g,oj) for 
small u, and in this case we can evaluate this quantity 
in the limit lo — > 0, 



x {[GT(P + ^e)GT{P,e)]-c.c) 
= TJ2 I dl[G^{p + q ,ie n )G^{p,\e n )] 



+ 
+ 



/"r(W-) 



(M x + My 

f " r{ ^'\vl+vl) 
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Q 2 } + 0(q 4 



XAl) 



where we have used -^-G^(p, \e n ) = [G^°\p, ie n )] 2 and 

TEsJ G ° ) (P^n)} n+lP = ^eg-/r(^), and Mj(j + 
x,y) is given below Eq. (42). 

On the other hand, contour C2 gives the imaginary part 
of Xa T q R (Q} w ) f° r sma U w <C Wf|q| as 



J 4tti 7 p p \ i/t -v p qj 
x i[Gfb,^ C ,.], + 0(, 2 ) 



p I 2r + i« p • g 



xiG^(p,k„) + 0( W 2 ) 



9 Jv 



-f'riZp,*-) 
2T/\q\ + iv p ■ q 



where we have used ^G^ R ^ A (p, e) 
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